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ABSTRACT 

The WMAP mission has produced sky maps from 7 years of observations at L2. We present the 
angular power spectra derived from the 7-year maps and discuss the cosmological conclusions that 
can be inferred from WMAP data alone. 

With the 7-year data, the temperature (TT) spectrum measurement has a signal-to-noisc ratio 
per multipole that exceeds unity for I < 919; and in band-powers of width Al ~ 10, the signal-to- 
noise ratio exceeds unity up to / = 1060. The third acoustic peak in the TT spectrum is now well 
measured by WMAP. In the context of a flat ACDM model, this improvement allows us to place 
tighter constraints on the matter density from WMAP data alone, ilmh^ 
epoch of matter-radiation equality, z, 

detected in the 7-year data with a significance of 20a, compared to 13a with the 5-year data. We now 
detect the second dip in the TE spectrum near I ^ 450 with high confidence. The TB and EB spectra 
remain consistent with zero, thus demonstrating low systematic errors and foreground residuals in 
the data. The low-Z EE spectrum, a measure of the optical depth due to reionization, is detected at 



0.1334t!] °[;^6^ and on the 
cq — 3196^]^'^. The temperature-polarization (TE) spectrum is 



5.5cr significance when averaged over I = 2-7: l{l + l)Cf ^/(27r) = 0.074 



1+0.034 
'0.025 



fjK^ (68% CL). We 



now detect the high-/, 24 < Z < 800, EE spectrum at over 8 a. The BB spectrum, an important 
probe of gravitational waves from infiation, remains consistent with zero; when averaged over I = 2-7, 
1(1 + 1)CP^ /{2tt) < 0.055 /xK^ (95% CL). The upper limit on tensor modes from p olarization data 
alone is a factor of 2 lower with the 7-year data than it was using the 5-year data (jKomatsu et al.l 
[20T0I ). 

The data remain consistent with the simple ACDM model: the best-fit TT spectrum has an effective 
of 1227 for 1170 degrees of freedom, with a probability to exceed of 9.6%. The allowable volume 
in the 6-dimensional space of ACDM parameters has been reduced by a factor of 1.5 relative to the 
5-year volume, while the ACDM model that allows for tensor modes and a running scalar spectral 
index has a factor of 3 lower volume when fit to the 7-year data. We test the parameter recovery 
process for bias and find that the scalar spectral index, n^, is biased high, but only by 0.09 a, while 
the remaining parameters are biased by < 0.15 a. 

The improvement in the third peak measurement leads to tighter lower limits from WMAP on the 
number of relativistic degrees of freedom (e.g., neutrinos) in the early universe: A^off > 2.7 (95% CL). 
Also, using WMAP data alone, the primordial helium mass fraction is found to be Fho = 0-28io i5, and 
with data from higher- r esolution CMB exper iments included, we now establish the existence of pre- 
stellar helium at > 3a (jKomatsu et al.l [20101 ). These new WMAP measurements provide important 
tests of Big Bang cosmology. 
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1. INTRODUCTION 

The Wilkinson Micr owave Anisotropy Probe (WMAP 
iBennett et all l2003bl lal) is a Medium- Class Explorer 
(MIDEX) satellite aimed at understanding cosmology 
through full-sky observations of the cosmic microwave 
background (CMB). The WMAP full-sky maps of the 
temperature and polarization anisotropy in five fre- 
quency bands provide our most accurate view to date of 
conditions in the early universe. The WMAP instrument 
is composed of 10 differencing asse mblies (DAs) spanniii g 
5 frequencies from 23 to 94 GHz (jBennett et al.ll2003b[ ): 
one DA each at 23 GHz (Kl) and 33 GHz (Kal), two each 
at 41 GHz (Q1,Q2) and 61 GHz (V1,V2), and four at 94 
GHz (W1-W4). Each DA is formed from two differen- 
tial radiometers which are sensitive to orthogonal linear 
polarization modes; WMAP measures both temperature 
and polarization at each frequency. The multi-frequency 
data facilitate the separation of the CMB signal from 
foreground emission arising both from our Galaxy and 
from extragalactic sources. The CMB angular power 
spectrum derived from these maps exhibits a highly co- 
herent acoustic peak structure which makes it possible 
to extract a wealth of information about the composi- 
tion and history of the universe, as well as the processes 
that seeded the fluctuations. 

With accurate measurements of the first few peaks 
in the angular power spectrum, CMB data have en- 
abled th e following advances in our understanding of cos- 
moloev (iSpergel et all [20031 [2007t iDunklev et all [200l 
iKomatsu et al.ll2009D : the dark matter must be non- 
baryonic and interact only weakly with atoms and ra- 
diation; the density of atoms in the universe is known 
to 3%, and accords well with Big Bang Nucleosynthesis; 
the measured acoustic scale at z = 1090, combined with 
the local distance scale and Baryon Acoustic Oscillation 
(BAG) data, demonstrates that the universe is spatially 
flat, to within 1%; the Hubble constant is determined to 
3% using only acoustic fluctuation data (CMB+BAO), 
and it accords well with local measurements; the pri- 
mordial fluctuations are adiabatic and Gaussian, and the 
spectrum is slightly tilted. 

The statistical properties of the CMB fluctuations 
measured by WMAP are close to Gaussian with ran- 
dom phase (IKomatsu et all 120031 : ISpergel et aD 120071 : 
IKomatsu et al.ll2009l ). There are several h ints of possible 
deviat i ons from this ca s e as d iscussed in IBennett et al.l 
(|2010[) : IKomatsu et"aI1 (|2010[) . If the fluctuations are 
Gaussian and random phase, then their statistical infor- 
mation content is completely determined by the angular 
power spectra of the sky maps. 

This paper derives the angular power spectra from the 
WMAP 7-year sky maps and presents the cosmological 
parameters that can be determined from them. The new 
results improve upon previous results in many ways: ad- 
ditional data reduces the random noise, which is espe- 
cially important for studying the temperature signal on 
small angular scales and the polarization signal on large 
angular scales; W band data is now incorporated in the 
TE spectrum measurement to improve precision; new 
simulations have been carried out to test the accuracy 
of parameter recovery and to test a model's goodness of 
fit. The result is the most accurate full sky measurement 
to date of CMB anisotropy down to an angular scale of 



- 0.25°. 

This paper is one o f six that acc o mpan y the 7-year 
WMAP data release. iJarosik et all (I2010D discuss the 
7-year map makin g process, s ystem atic error limits, 



and basic results 
foreground emission . 



Gold et all ()2010l ) discuss galactic 
and its removal in 7-year data. 



IBennett et al.l (|2010f ) discuss possible anomalies in the 
WMAP CMB maps. IKomatsu et"aI1 (|2010[) discuss the 
interpretation of the WMAP data , in combinatio i i with 
other relevant cosmological data. IWeiland et al.l ()2010l ) 
discuss the WMAP measurements of the outer planets 
and selected bright sources for use as microwave calibra- 
tors. 

The layout of this paper is as follows. In Sjll we present 
the WMAP 7-year power spectra. In Jj3l we discuss sim- 
ulations that were performed to test for bias in our cos- 
mological parameter fits. In Sj4l we discuss cosmological 
conclusions that can be drawn from WMAP data alone, 
and in SjSjwe discuss the goodness of fit of the 6-parameter 
ACDM theory. We conclude in 

2. SEVEN-YEAR POWER SPECTRA 

In this section we present the temperature and polar- 
ization power spectra derived from the 7- year sky maps 
and compare them to the 5-year spectra. 

2.1. Definitions and Methodology 

Since WMAP measures both temperature and polar- 
ization, there are multiple power spectra to consider. On 
a sphere, the temperature field can be decomposed into 
spherical harmonics. 
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where n is a unit direction vector and ot refers 
specifically to the temperature field. Likewise, the 
Q and U Stokes parameters for linear polarization 
can be de compose d into complex spin-2 ha rmonics 
(|Newman fc PenrosilT96l iGoldberg et al.lll967| ). 
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The spin-2 coefficients can then be combined to repre- 
sent polarization modes that have no curl (E modes) 
and modes that have no diver gence (B modes). These 
are given by the coefficients (jZaldarriaga fc Seliaklll997l : 
lLarsonll200l . 
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(jKamionkowski et al.l ()1997[ ) use an alternative ap- 
proach.) The angular power spectra are related to these 
modes according to 
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where X, Y = T, E, or B. The data are cur- 
rently consistent with being isotropic and Gaussian dis- 
tributed, but this co ndition should continue to be tested 
([Bennett et alJl2010f ). 

There are 6 independent power spectra than can be 
constructed from the temperature and polarization data, 
TT, TE, TB, EE, EB, and BB, though in theories in 
which parity is conserved, TB and EB are expected to 
be zero. In general, foreground signals (and systematic 
effects) can produce non-zero TB, EB, and BB so these 
spectra provide a good test for residual polarization con- 
tamination. 

For the 7-year analysis, we use the same combina- 
tion of estimators that were used with the 5-year data 
(|Nolta et al.ll2009| ). This combination is a trade-off be- 
tween statistical accuracy and computational speed. We 
present new tests of the accuracy of the likelihood func- 
tion constructed from these estimators in ^ To summa- 
rize the combination: for low-/ TT (< 32) we compute 
the likelihood of a model directly fro m the 7- year Inter - 
nal Linear Combination (ILC) maps (IGold etal.1 (20101) . 
For high-/ (/ > 32) TT we use the MA STER pseudo-C; 
quadratic estimator (jHivon et al.ll200^ . For low-/ polar- 
ization, / < 23 TE, E E, fc BB, we use the pixel-space 
estimator described in iPage et ahl ()2Q07t) , and for high-/ 
TE (/ > 23) we use the MASTER quadratic estimator. 

2.2. Changes Affecting the 7-Year Spectra 

Several data processing and analysis changes were ap- 
plied to the 7- year data which resulted in improvements 
beyond those which would be expected from additional 
integration time. 

2.2.1. Map-Making with Asymmetric Masking 

A new map-making technique was adopted for the 7- 
year data which combines o ptimal noise handlin g with 
"asymmetric" data masking (jJarosik et al.l[20To| ). With 
this change, certain regions in the 7-year maps employ 
more data samples than they would have with the pre- 
vious pipeline. These "Galactic echo" regions are thus 
more sensitive than a simple 5-year to 7- year integration 
time scaling would predict. 

2.2.2. Multipole Range 

The additional sensitivity afforded by more data has 
made it possible to extend the usable multipole range 
in the power spectra and likelihood code. For the TT 
data we extend the upper multipole limit, /max, from 
1000 to 1200. For the TE spectrum, we have determined 
that high-/ W-band polarization data are sufficiently free 
from systematic effects that they can be em ployed in the 
TE spectrum estimate (jJarosik et al.|[20lol ). This signifi- 
cantly improves the sensitivity in the 7- year TE spectrum 
and allows us to extend the TE multipole limit from 450 
to 800. 

2.2.3. Mask and f^ky 

The 7-year sky masks have been augmented slightly 
using a analysis of the Q-V & V-W difference maps, 
after the normal template cleaning had been applied 
(jJarosik et al.ll2010l: iGold et al.ll2010D . This results in a 
slightly more conservative mask which decreases the un- 
masked sky fraction by ~ 3% (from 81.7% to 78.3% for 



the KQ85 cut — the new cut is denoted KQ85y7). Given 
the threshold applied during the construction of the 
extended mask, residual foreground signals outside the 
mask are essentially undetectable on the scale of the in- 
strument noise in a ^2° pixel, so the data are fractionally 
more robust to foreground contamination. 

The power spectrum sensitivity depends on sky cut 
according to AC/ oc f~^^ where /sky is (approximately) 
the fr action of sky that survives the cut (jHinshaw et al.l 
120031 IVerde et alJl2003| ). In practice, /sky is a function 
of / that is calibrated with simulations, so there is a dif- 
ferent constant of proportionality at each /, but it scales 
with the fraction of usable sky area. Thus the increased 
sky mask results in a slight loss of sensitivity in the TT 
spectrum for / < 550, where the spectrum is sky variance 
limited. 

For the TE spectrum, we have generated new simula- 
tions to test the calibration of /sky.TE which enters into 
the error propagation from sky maps to spectra. The 
mean deduced from our simulations was 760 for spec- 
tra with 777 degrees of freedom, a factor of 1.022 too low, 
indicating that our previous TE error estimate was a fac- 
tor of 1.011 too high. Therefore we have scaled /sky,TE 
by 1.011 to produce a unit mean psr degree of free- 
dom, following th e precedent used to calibrate /sky.TT 
from simulations (jVerde et al.l |2003() . In ^J5] we present 
the TE of the 7-year flight data and conclude that 
the ACDM model fits the TE data weU. 

2.3. Temperature (TT) Spectrum 

For / < 32, the spectrum is obtained using a Blackwell- 
Rao estimator appl i ed to a chain of Gibbs samples 
iWandeh et al."2004|; iJewell et al.ll200l IChu et al.ll2005l: 
[Dunklcy et al. 200!|) based on the 7-year ILC map and 
the KQ85y7 mask. T he specifications used to sample the 
map are described in lDunklev et al.l (j2009l ). For / > 32, 
the spectrum was derived from the MASTER pseudo- 
Ci quadratic estimator applied to the 7- year, te mplate- 
cleaned V and W-band maps (jCold et al.ll2010l ). (The 
MASTER spectrum is technically derived from / = 2 
to 1200, then the / = 2 — 32 portion is discarded, but 
correlations induced by mode coupling are retained for 
/ > 32.) The pseudo-a;,„ coefficients were computed from 
the A^sidc = 1024 maps for each single year and each 
single-DA, V1-W4. For / < 600 the coefficients were 
evaluated with uniform pixel weights, while inverse-noise 
weights were used for / > 600. (The transition was made 
at / = 500 in the 5-year analysis.) As noted above, we 
adopt a slightly larger sky mask, denoted KQ85y7. 

The pseudo-C/ cross-power spectra are computed from 
all off-diagonal pairs of pseudo-a/m coefficients, 

' ^ 2/ + 1 ^ ^^im^Im IDJ 
m— — / 

wher e i, j refer to a DA- year combination (jHinshaw et al.l 
I2007f ) and the tilde indicates a pseudo-quantity (spec- 
trum or coefficient). These component pseudo-spectra 
are deconvolved using the MASTER formafisnH, and 
the results are combined by band into VV, VW, and 

In principle, we can obtain a modestly more sensitive 
spectrum estimate at intermediate multipoles by employing 
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Figure 1. The 7- year temperature (TT) power spectrum from WMAP. The third acoustic peak and the onset of the Silk damping tail 
arc now well measured by WMAP. The curve is the ACDM model best fit to the 7-year WMAP data: n,,h^= 0.02270, Vlch^= 0.1107, 
f2A= 0.738, T= 0.086, ns= 0.969, A|j= 2.38 x 10^^, and j4g2= 0.52. The plotted errors include instrument noise, but not the small, 
correlated contribution due to beam and point source subtraction uncertainty. The gray band represents cosmic variance. A complete error 
treatment is incorporated in the WMAP likelihood code. The points are binned in progressively larger multipole bins with increasing I; 
the bin ranges are included in the 7-year data release. 

spectra evaluated with the KQ85y7 mask. (Most of the 
cosmological parameters reported in this paper were fit 
using a prehminary version of the hkehhood that had a 
smaU masking error that produced a shghtly biased TT 
spectrum at high-Z and a correspondingly higher resid- 
ual source amplitude, which mostly compensated for the 
bias. We have checked that substituting the correct TT 
spectrum has a negligible effect on the parameter fits.) 
After this source model is subtracted from each band, 
the spectra are combined to form our best estimate of 
the CMB signal, shown in Figure [T] 
1100 1200 The 7-ycar power spectrum is cosmic variance limited, 
i.e., cosmic variance exceeds the instrument noise, up to 
I = 548. (This limit is slightly model dependent and can 
vary by a few multipoles.) The spectrum has a signal- 
to-noise ratio greater than one per ^-mode up to Z = 919, 
and in band-powers of width A/ = 10, the signal-to- noise 
ratio exceeds unity up to I = 1060. The largest improve- 
ment in the 7-year spectrum occurs at multipoles I > 600 
where the uncertainty is still dominated by instrument 
noise. The instrument noise level in the 7-year spectrum 
is 39% smaller than with the 5-year data, which makes it 
worthwhile to extend the WMAP spectrum estimate up 
to / = 1200 for the first time. See Figure[2]for a compari- 
son of the 7- year error bars to the 5-year error bars. The 
third acoustic peak is now well measured and the onset 
of the Silk damping tail is also clearly seen by WMAP. 
As we show in i|4l this leads to a better measurement 
of flmh^ and the epoch of matter-radiation equality, Zcq, 
which, in turn, leads to better constraints on the effective 
number of relativistic species, A^off, and on the primor- 
dial helium abundance, Ihc- The improved sensitivity 



Figure 2. The high-i TT spectrum measured by WMAP, showing 
the improvement with 7 years of data. The points with errors use 
the full data set while the boxes show the 5-year results with the 
same binning. The TT measurement is improved by >30% in the 
vicinity of the third acoustic peak (at / 800), while the 2 bins 
from I = 1000-1200 are new with the 7-year data analysis. 

WW spectra for the purposes of removing the residual 
point source amplitude. The unresolved point source 
contribution to the sky continues to be treated as a 
power law in thermody namic te mperature, falling as 
(|Nolta et al.ll20Q9D . but see iColombo fc Pierpaolil 
(|2010D for an alternative approach to the spectral de- 
pendence. Using the same fitting methodology as in the 
5-year analysis, we find its amplitude to be 10'^j4ps = 
9.0±0.7//K2 sr, when fit to the 7-year Q, V, and W band 

weighting in the computation of the pseudo-aj,„. We are currently 
developing such code for use in cross-power spectra with the inten- 
tion of applying it to the final 9-year data. 
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Figure 3. The 7-year temperature-polarization (TE) cross-power 
spectrum measured by WMAP. The second trough (TE<0) in the 
spectrum in the vicinity of I = 450 is now clearly detected. The 
green curve is the ACDM model best fit to the 7- year WMAP data, 
as in Figure [T] The plotted errors depict the diagonal elements 
of the covariance matrix and include both cosmic variance and 
instrument noise. A complete error treatment is incorporated in 
the WMAP likelihood code. Note that the plotted spectrum is 
(Z -I- l)C™/(27r), and not l{l + 1)C™/{2tt). 

at high I is also important for higher-resolution CMB 
experiments that use WMAP as a primary calibration 
source. 

2.4. Temperature- Polarization (TE, TB) Cross Spectra 

The 7- year temperature-polarization cross power spec- 
tra were forme d using the same methodology as t he 5- 
year spectrum ([Page et al.l 120071: iNolta et al.i(2009l ). For 
/ < 23 the cosmological model likelihood is estimated 
directly from low-resolution temperature and polariza- 
tion maps. The temperature input is a template-cleaned, 
co-added V-l-W band map, while the polarization in- 
ut is a ternplate- cleaned, co-added Ka-|-Q-|-V band map 
Gold et al.l 120091 ) ■ In this regime, the spectrum can be 



inferred from the conditional likelihood of Ci values (in- 
dividual or binned), but these estimates are only used 
for visualization. 

For ^ > 23, the temperature-polarization spectra are 
derived using the MASTER quadra tic estimator, ex - 
tended to include polarization data (jPage et al.l I2007D . 
(As above, the MASTER spectrum is evaluated from 
I = 2, but the result from Z = 2 — 23 is discarded.) The 
temperature input is a template-cleaned, co-added V-l-W 
band map, while the polarization input is a template- 
cleaned, co-added Q-l-V+W band map. The inclusion of 
W-band data in the high-/ TE and TB spectra is new 
with the 7-ycar data release (jJarosik et al.|[20Tol ). Since 
the W band radiometers have the highest angular resolu- 
tion, the inclusion of W band significantly enhances the 
sensitivity of these high-/ spectra. 

The 7- year TE spectrum measured by WMAP is shown 
in Figure [31 For all except the first bin, the MAS- 
TER values and their Gaussian errors are plotted. The 
first bin shows the conditional maximum likelihood value 
based on the pixel likelihood mentioned above. The 
slight adjustment for /sky,TE is included in the error bars. 
With two additional years of integration and the inclu- 
sion of W band data, we now detect the TE signal with 
a significance of 20cr, up from 13(t with the 5-year data. 
Indeed, for 10 < / < 300, the TE error is less than 65% of 



the 5-year value, and for I > 300 the sensitivity improve- 
ment is even larger due to W band's finer resolution. At 
/ = 800 the 7-year TE error is 36% of the 5-year value. 
A qualitatively new feature seen in the 7-year spectrum 
is a second trough (TE<0) near I = 450. See Figure |4] 
for a comparison of the 7-year to 5-year error bars, for 
the TE and TB spectra. Overall, the TE data are quite 
consistent with the simplest 6-parameter ACDM model; 
we discuss its goodness-of-fit in ^ 

The observed TE signal is the result of a specific polar- 
ization pattern around hot and cold spots in the temper- 
ature anisotropy. In particular, the acoustic peak struc- 
ture in TE corresponds to a series of concentric rings of 
alternating radial and tangen tial polarization (relat ive to 
a radial reference direction). iKomatsu et al.l ( 20101) per- 
form a stacking analysis of the 7-year temperature and 
polarization maps and show that the effect is detected in 
the 7-year WMAP sky maps with a significance of 8cr. 

The 7-year TB spectrum measured by WMAP is shown 
in Figure [5l In this case, because the signal-to-noise ra- 
tio is low, the MASTER points and their Gaussian errors 
are plotted over the full I range, including the first bin. 
The measured spectrum is consistent with zero: the 
for the null hypothesis (TB=0) is 793.5 for 777 degrees 
of freedom. The probability to exceed that amount is 
33%. The absence of a detectable signal is consistent 
with the ACDM model, which predicts zero. It is also 
an indication that systematic errors and foreground con- 
tamination are not significant at the level of ^ 0.1 /xK^ 
i n {I + DCP 



IKomatsu et~all (|2010[) use the 7-year TE and TB data 
to place limits on polarization rotation due to parity vio- 
lating effects. Polarization rotation would cause TE sig- 
nal generated at the last scattering surface to transform 
into observed TB power. The absence of TB signal leads 
to an upper limit on rotation of Aa = — l.l°±1.4°(stat)± 
1.5°(sys). 

2.5. Polarization (EE, EB, BB) Spectra 

We begin by discussing the low-Z polarization spectra, 
and then move on to the high-Z EE spectrum. 

The most reliable way to estimate the low-/ polar- 
ization spectra is to use the pixel-space likelihood code 
to generate the posterior distributions of individual (or 
binned) C/ values. In the 7- year data, this code is based 
on a co-added Ka-|-Q-|-V map. The most conservative, 
but costly, method is to produce a Markov Chain that 
allows each C; to vary independently; the resulting dis- 
tribution of any single C/ will be the marginalized dis- 
tribution for that multipole moment. A Gibbs sampling 
technique could also be used, but this works best with 
a high signal-to-noisc ratio. However, Gibbs sampling in 
lower signal-to-nois e regions can be p erformed success- 
fully, as shown by iJewell et al.l ()2009l ) . A much more 
tractable approach is to compute the conditional likeli- 
hood in which the likelihood of a single Ci is evaluated 
while all other moments are held fixed. We adopt the 
latter approach to visualize the low-Z EE and BB spec- 
tra. This method has also been used in previous WMAP 
papers as well as (for example) iGruppuso et al.l ()2009l ) 
in their verification of the 5-year WMAP low-/ spectra. 
In the context of parameter fitting, the estimated C; are 
constrained to vary according to the model. 

Figure [6] shows the conditional likelihood for the EE 
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Figure 4. The TE and TB high-( spectra measured by WMAP, showing the improvement with 7 years of data. The points with errors 
use the full data set while the boxes show the 5-year results with the same binning. The spectra are greatly improved by the addition of 
W-band data. The non-detection of TB signal is expected; it provides a good che ck of sys tematic errors and fore ground residuals, and can 
be also used to set limits on polarization rotation due to parity- violating effects ( i|2.4l and lKomatsu et aLl 1120101 )). 

the C]^^ are fixed to the best- fit ACDM model for I' ^ I. 
The red curves are tlie analogous distributions when the 
reference spectrum is taken to be the maximum likeli- 
hood spectrum. This maximum likelihood spectrum was 
obtained by numerical maximization of the likelihood 
code for the TT, TE, EE, and BB spectra for 2 < Z < 10, 
a maximization in 36 dimensions, while the spectra at 
/ > 10 were fixed at the best-fit ACDM model. Save 
for I — 3 and 6, the likelihood curves are relatively in- 
sensitive to the difference between these two reference 
spectra. From these curves it is clear that the majority 
of the statistical weight in the \ow-l EE detection is at 
I = 4, with / = 2 also contributing significant power. 
A standard reionization scenario would give rise to a 
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Figure 5. The 7- year temperature-polarization (TB) cross-power 
spectrum measured by WMAP. This spectrum is predicted to be 
zero in the basic ACDM model and the measured spectrum is con- 
sistent with zero. TB provides a useful null test for systematic 
errors and foreground residuals. IKomatsu et al.l I I2010I) use the TB 
and TE spectra to place an upper limit on polarization rotation 
due to parity- violating effects. The TB fo'' the null hypothe- 
sis (TB=0) is 793.5 for 777 degrees of freedom. The probability 
to exceed that amount is 33%. Note that the plotted spectrum is 
(Z -1- l)C™/(27r), and not l{l + l)C™/{2n). 

multipoles from I = 2-7 for two different reference spec- 
tra. The black curves show the likelihood of C?^ when 



relatively flat spectrum in Cf^ = l{l + 1)0^^/(2^) over 
the range Z = 2 — 7, so it is of interest to evaluate the 
posterior distribution of a band power with constant C 
over this range. As shown in Figure [71 we find 



C™7 = 0.074 



+°:^^4 (68% CL). 



(7) 



This result was obtained with the pixel likelihood code, 
and so the error bars include cosmic variance. Addition- 
ally, a model with zero TE and EE power for / ~ 2-7 
is disfavored at 5.5a relative to the most-likely constant 
band-power in this / range. 

Figure [8] shows the conditional likelihood for the BB 
multipoles from I = 2-7 for two different reference spec- 
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Figure 6. Conditional likelihoods of l{l + l)Cp^/(27r) for / = 2-7, computed with the WMAP UkeUhood code. The Ka, Q, and V bands 
contribute to the \ow-l polarized pixel likelihood. In each panel, the black curve is the conditional likelihood for a given multipole when all 
other multipoles are held fixed at the value of the best-fit ACDM model (indicated by the black diamonds). The red curve in each panel is 
the conditional likelihood when all other multipoles are held fixed at the maximum likelihood value of the spectrum, indicated by the red 
diamonds. The maximum likelihood spectra were determined by a numerical maximization of the WMAP likelihood code, for I = 2-10, 
for TT, TE, EE, and BB. Points with I > 10 were fixed at the best-fit ACDM value. 



tra. The black curves show the likchhood of when 
the C;?^ are fixed to the best fit ACDM model (zero, ex- 
cept for small contributions from lensing) for I' I. The 
red curves are the analogous distributions when the ref- 
erence spectrum is taken to be the maximum likelihood 
spectrum for I < 10 and the ACDM spectrum (again, 
effectively zero) for ^ > 10. Save for / = 2, the pos- 
terior likelihood curves arc insensitive to the reference 
spectrum. There is no significant detection of BB power 
in any single multipole in the 7-year WMAP data. As 
shown in Figure [HI we evaluate the posterior likelihood 
of a single constant band-power from 1 = 2 — 7 and find 
it is also consistent with zero. We place an upper limit 
of 

Cf^T < 0.055 fiK^ (95% CL) (8) 
using the 7-year WMAP data, which is more than a fac- 



tor of 2 lower than the 5-year limit of 0.15 /iK^. 

The high-Z EE spectrum is constructed using the po- 
larized spectra that were used for the TE spectrum, dis- 
cussed in ^ 32.41 FigurcfTOlplots the 7- year WMAP data on 
top of the ACDM EE spectrum that best fits the WMAP 
data. Note that the high-Z EE spectrum is not included 
in the likelihood code, so the theory curve is not a best 
fit to the high-Z EE spectrum. 

Using error bars with cosmic variance derived from the 
best fit ACDM model, we find = 830.6, for the 777 
degrees of freedom in the multipole range 24 < I < 800. 
The probability to exceed this value is 8.9%, which 
is low, but not significantly so. For a model with no EE 
spectrum, = 897.3. The difference is Ax^ = 66.7, 
which is just over an 8 cr detection of the high-Z EE spec- 
trum. 

The data prefer an EE spectrum with higher ampli- 
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Figure 7. The relative for a single bin of l{l + l)Cp^/{27r) 
from I = 2-7, conditioned on the best-fit ACDM spectrum. For 
simplicity, wc set Cj^ = for I = 2-7 for this evaluation, so that 

the constraint < '■JcJ'^Cf^ will always be satisfied. The 

vertical red lines indicate where Ax^ = 1, 4, 9, and 16, corre- 
sponding to 1, 2, 3, and 4 a confidence limits on EE. The EE=0 
point in this bin has Ax^ = 26.5; additionally, setting TE=0 in the 
2-7 bin raises by 3.5 relative to the best-fit ACDM TE spec- 
trum. Thus the full change in between the best-fit model and 
the EE=TE=0 model is 30, corresponding to a 5.5cr detection of 
EE power in this bin, with 1(1 + l)Cf '=/{27r) = OmAtoUt fiK^ 
(68% CL). 

tudc than the best-fit ACDM model. To quantify this, 
we find the seale faetor a that causes the theory EE 
spectrum to best fit the data, over the multipole range 
24 < Z < 800. The scale factor is a = 1.52 ± 0.21. This 
indicates a 2.5 a preference for a higher amplitude EE 
spectrum, which we consider to be worth further inves- 
tigation, but do not believe to be a significant deviation 
from the ACD M theory. Note tha t other experiments, 
such as QUaD (jBrown et al.ll2009[) . find the EE spec- 
trum to be consistent with the prediction of the best-fit 
model. 

To verify that we are not seeing the power spectrum 
in just one band and not the others, we take difference 
spectra among the Q, V, and W bands. Figure [TT] shows 
two sets of difference spectra. These spectra are consis- 
tent with zero, as expected, and demonstrate that the EE 
power spectrum is present in all three frequency bands. 

2.6. WMAP Likelihood Code 

Before discussing cosmological parameter fits in the re- 
mainder of the paper, we review the WMAP likelihood 
code which forms the basis for the fits. 

The basic structure of the likelihood code is unchanged 
from the 5-year release. For the / > 32 temperature 
data, the model spectrum is compared to the MASTER 
spectrum, described above, using a Gaussian plus log- 
normal app r oxima tion to the likelih o od, a s described in 
iBond et ail ([l998l ) and iVerde et"all ([200l . For I < 32 
a Blackwell-Rao estimator is used to determine the like - 
lihood of a model TT spectrum (jDunklev et al.l 120091 ). 
This estimator encodes both the low-/ spectrum and an 
accurate description of its no n- Gaussian errors. It is con- 
structed from a set of Gibbs samples that contain power 
spectra and C MB maps that a re statisti cally consistent 
with the data (|Wandelt et al.l [2004; .Jewell et al.l 120041: 



lEriksen et aIll2004D . The 7-year input to the Gibbs chain 
mimics the 5-year input: we smooth the Internal Linear 
Combination (ILC) map to 5° (Gaussian FWHM); de- 
grade it to resolution A'sido = 32; and add Gaussian white 
noise with 2 fiK rms to each pixel. The data are masked 
with the KQ85y7 mask degraded to A'sido = 32, and the 
Gibbs sampler is run to produce input for the Blackwell- 
Rao estimator; further details are given in lDunklev et all 
(2009). 

For the polarization data, we use a similar hybrid 
scheme: for the / > 23 TE data we compare a model 
spectrum to the MASTER spectrum using a Gaussian 
likelihood. (TE is the only high-/ polarization data used 
in the WMAP likelihood code.) For / < 23, the likelihood 
of model TE, EE, and BB spectra are obtained using a 
pixel-space likelihood which is based on the A"sidc = 8 
map mentioned in N2.4I and described in ([Page et al.l 
l2007h . 

The likelihood code includes several important fac- 
tors: mode coupling due to sky masking and non-uniform 
pixel weighting (due to non- uniform noise); beam win- 
dow function uncertainty, which is correlated across the 
entire spectrum; and residual point source subtraction 
uncertainty, which is also highly correlated. The treat- 
ment o£_these_effectsis ur^^ 5-year anal- 
ysis (|Nolta et al.ll2009l: iDunklev et al.ll2009D . 

Note added in revision — The results in this paper were 
prepared using version 4.0 of the WMAP likelihood. 
Since the initial submission of this paper, two small errors 
in the likelihood code came to light. 1) The original com- 
putation of the TT spectrum used an incorrect monopole 
subtraction which resulted in a small amount of excess 
power at high /, and a corresponding elevation of the 
best-fit residual point source amplitude. Correcting the 
monopole subtraction reduced the high-/ power slightly 
which produced a correspondingly lower residual point 
source amplitude, from 11.0 x 10""^ to 9.0 x 10~^^K^ sr. 
2) Due to a simulation configuration error, the TE /sky 
recalibration factor used in version 4.0 was 3.8% larger 
than the final value reported in i i2.2.3l The first of these 
changes will not affect the simulations in fJ31 because 
they lack a monopole, and the new value for TE /sky 
has been used for the ^parameter recovery simulations. 
The goodness of fit statistics for the TT and TE spectra 
in Sj5] compare the best fit ACDM theory spectrum from 
the version 4.0 Markov chains (with RECFAST version 
1.4.2) to the version 4.1 likelihood data. See Appendix 
[B]for a comparison of parameters, when estimated with 
the original and updated versions of the code. 

2.6.1. CAMB 

For computing theoretical power spectra, we use the 
Code for Anisotropics in the Microwave Background 
(CAMB) (|Lewis et aLll2000[) which is b ased on the earlier 
code CMBFAST (ISeliak fc Zaldarriaea ,1996.l FI Since 
2008 the CAMB package has supported improved mod- 
eling of reionization, as follows: 1) inclusion of helium 
reionization, assuming that helium is singly reionized 

^"^ Wc use the November 2008 version of CAMB, which was up- 
dated to remove a bug affecting lensed non-flat models in February 
2009. The update was inadvertently not included in our current 
analysis, but we have checked that the effect on spectra is at the 
sub-percent level. 
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Figure 8. Conditional likeliiioods of l{l + l)Cf^ /{2-k) for I = 2-7, computed witli the WMAP likcliliood code. Tlic Ka, Q, and V bands 
contribute to tiie low-Z polarized pixel likelihood. In each panel, the black curve is the conditional likelihood for a given multipole when 
all other multipolcs are held fixed at the value of the best-fit ACDM model (effectively zero, except for gravitational lensing) as indicated 
by the black diamonds. The red curve in each panel is the conditional likelihood when all other multipolcs are held fixed at the maximum 

The maximum likelihood spectrum was determined as stated in the 

cal parameters. The fractional lowering is largest at high 
/ and is about 1% at / = 1000. 

We discuss the effect of these changes in Appendix [b1 
The parameter result s presented in this paper, and in 
iKomatsu et"all (|2010[ ). use the original RECFAST ver- 
sion 1.4.2. 

3. PARAMETER RECOVERY BIAS TESTS 

In this section, we describe a test for bias in the WMAP 
parameter recovery process. The parameters of the basic 
ACDM model are: the physical baryon density, Vl^h?] 
the physical cold dark matter density, ftch^; the dark 
energy density, in units of the critical density, the 
amplitude of primordial scalar curvature perturbations 
at fc = 0.002 Mpc~^, A^; the power-law spectral index 
of primordial density (scalar) perturbations, n^; and the 
reionization optical depth, r. The above parameters are 



likelihood value of the spectrum, indicated by the red diamonds, 
caption to Figure [6] 

at the same time as hydrogen and doubly reionized at 
z ~ 3.5 (this slightly lowers the redshift of reionization 
for a given optical depth); 2) the width of reionization 
can be varied without changing the optical depth. We 
use a width Az = 0.5 as standard. Seven-year reioniza- 
tion results are discussed in ^4.2.71 

Shortly after the 7-year WMAP data were released, a 
new version of CAMB was made available, incorporating 
an updated version of the code used to m odel r ecombi- 
natio n: RECFAST (ISeager et al.l[T999l l2000: Won g et al.l 
120081 iScott fc Mossll2009D . The parameter recovery was 
run with RECFAST version 1.4.2, and the January 2010 
version of CAMB updates this to RECFAST version 1.5. 
The primary change is in the optical depth, due to more 
accurate modeling of the physics of the hydrogen atom. 
Figure [12] shows how much the new version of RECFAST 
lowers the power spectrum, for a given set of cosmologi- 
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Table 1 

Cosmological Parameter Definitions^ 



Parameter 



Description 



Fit parameters 



As2 



Physical baryon density 
Physical cold dark matter density 

Dark energy density {w = —1 unless otherwise noted - see below) 
Amplitude of curvature perturbations, fco = 0.002 Mpc^^ 
Spectral index of density perturbations, fco = 0.002 Mpc^^ 
Reionization optical depth 

Amplitude of the Sunyaev-Zel'dovich spectrum'' 



Derived parameters 



to 
Ho 



Age of the universe (Gyr) 

Hubble parameter. Ho = lOO/i km s"^ Mpc"'^ 

Amplitude of density fluctuations in linear theory, 8 Mpc scale 
Redshift of matter-radiation equality 
Redshift of reionization 



Extended parameters 

drii, / d\nk Running of scalar spectral index 

T Ratio of tensor to scalar perturbation amplitude, ko = 0.002 Mpc"'^ 

«_i Fraction of anti-correlated CDM isocurvature modes (se e i|4.2.3ll 

ao Fraction of uncorrelated CDM isocurvature modes (see i|4.2.3p 

fife Spatial curvature, Sl^ = 1 — fitot 

w Dark energy equation of state, w = Pde/pde 

N^fi Effective number of relativistic species (e.g., neutrinos) 

^Ho Primordial Helium fraction, by mass 

Az Width of reionization (new parameter in CAMB, see i|4.2.7|l 

^ Cosmological parameters discussed in this paper. A complete tabulation of the 
marginalized parameter values for each of the models discussed in this paper may be 
found at http://lambda.gsfc.nasa.gov 

^ The Sunyaev-Zel'dovich (SZ) amplitude is not sampled in the parameter recovery 
simulations, because the SZ effect is not included in the simulation power spectrum. 
Asz is sampled in the 7-year WMAP Markov c hain s, with a flat prior < Agz < 2, 
but is unconstrained by the WMAP data. See i|4.1l 



Binned BB, / -- 




0.00 0.05 0.10 0.15 

Z(Z+1)C,BB/27I [(XK2] 



0.20 



Figure 9. The relative x'^ for a single bin of Z(i-f l)Cp^/(27r) from 
I = 2 — 7, conditioned on the best-fit ACDM spectrum. The vertical 
red lines indicate where Ax^ = 1, 4, 9, and 16, corresponding to 
1, 2, 3, and 4 a confidence limits on BB. We find an upper limit of 
l{l + l)CpB/(27r) < 0.055 AtK2 (95% CL). 

sampled with flat priors and are sufficiently constrained 
by the WMAP data that boundaries to these priors do 
not have to be specified. Nevertheless, our Markov chain 
code adds the following constraints: 0.001 < Vlbh? < 0.2, 
0.0 < Qch^ < 0.5, < f^A, 0.04 < h"^ < 1.0, and 0.01 < 
T < 0.7. In this model, the Hubble constant, Hq = 100ft- 
km s~^ Mpc^"'^, is implicitly determined by the flatness 
constraint, Qb + flc + = 1- The Sunyaev-Zel'dovich 
(SZ) effect is not included in these simulations, nor is 
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Figure 10. WMAP detects the high-Z EE power spectrum. The 
green curve is the best fit ACDM model power spectrum, and the 
data are a combination of Q, V, and W band data. In the multipole 
range 24 < Z < 800, the detection is above 8(t. 

the Agz parameter sampled; see H4.ll for details on this 
parameter. Table [1] gives a description of the parameters 
considered in this paper including both fundamental and 
derived quantities. 
We generate 500 simulations of WMAP multi- 
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Figure 11. The Q-W and V-W band difference spectra for EE 
arc consistent with zero, as expected. The green line is the ACDM 
EE spectrum, which should not match this data, but is plotted 
to illustrate the size of the null-test error bars, compared to the 
detected signal. 
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Figure 12. The red line in this figure represents the percentage 
change in the C^-^-^ spectrum, due to the change in RECFAST 
from version 1.4.2 to version 1.5, including additional physics of 
the hydrogen atom. 

frequency sky map data with known cosmological pa- 
rameters, then verify that we recover the correct param- 
eters from those data. To our knowledge, this is the 
first statistical test of the likelihood in which multiple 
independent realizations are combined to test for bias at 
the ~ 0.1a level. Since we start the test with simulated 
sky maps, this work tests: the MASTER deconvolution 
of the masked pseudo-power-spectra; error propagation 
from maps to parameters; the code for combining V and 
W band data into a single TT spectrum; the code for 
combining Q, V, and W band data into a single TE spec- 
trum; the low-^ pixel-space likelihood codes for temper- 
ature and polarization; and the algorithm for combining 
these hybrid inputs into a single likelihood per model. 

There have been previous studies of the accura cy of 
the WMAP likelihood code. lO'Dwver etHI (|2004[ ) per- 
formed a Bayesian analysis of the 1-year WMAP data 
and found temperature p ower spectra la r gely c onsistent 
with those reported by iHinshaw et al.l (|2003| ). How- 
ever, they pointed out that the MASTER algo rithm does 
not ac curately represent the errors at low /. iChu et ahl 
()2005[) investigated cosmological parameters using the 
statistically exact BlackwcU-Rao estimator as part of 
the likelihood code at low I, and found shifts of up to 



0.5(7, compared to th e MASTER algorith m. These is- 
sues were addressed bv lSpergel et al.l feOOTi ) in the 3-year 
WMAP analy sis by using an A'sidc = 8 pixel likelihood. 
lEriksen et al.l ([2007) pointed out that the Ngidc = 8 code 
biased Cf^ slightly high in the range 12 < I < 30, 
which in turn biased n,^ slightly lo w. As a result, the 



final version of iSpergel et al.l (|2003| ) used an TVside = 16 
by code. Since then, the pixel likelihood code has been 
replaced with the Blackwell-Rao estimator, which accu- 
rately describes the C^"^ power spectrum up to / = 32 
(jPunklev et al.|[2009l ). 

The focus on arises because simple inflation mod- 
els predict its value to be slightly less than 1 (typically 
^0.96, which is termed spectral "tilt") and the best-fit 
value from previous WMAP analyses is in that range. 
However, the uncertainty is such that the deviation from 
1 is about 3cr, so small changes in the best-fit value can 
alter one's interpretation of significance. By comparing 
the WMAP lik elihood code to a Gau ssianized Blackwell- 
Rao estimator, iRudiord et al.l (|2009[ ) report a bias in Ug 
of -1-0. 6ct, which reduces the evidence for spectral tilt. In 
the results reported below, we do not find evidence for 
such a bias in the WMAP likelihood. In particular, using 
the 7-year WMAP data, we find = 0.963 ± 0.014 for 
a ACDM model fit. 

The pipeline for the WMAP data has already been ex- 
tensively tested by the WMAP team. This testing was 
in progress during the planning phase of the mission, and 
the pipeline continued to be refined after launch. Power 
spectrum reco nstruction from maps was simulated for the 
first year data ([Hinshaw et al.l2003() . The likelihood code 
was calibrated (by adjusting the effective sky fraction, 
/sky.Tx) with 100,000 simulations, so that the values 
reported from the likelihood could be us ed for goodness 
of fit tests as well as m odel comparison (jHinshaw et al.l 
I2OOI I Verde 'etsIMOO^ . Simulations of multiple years of 
time ordered data have shown that maps can be recon- 
structed, the WMAP data calibrated from the annual 
dipole modulation, and correct maps of the microwave 
sky recovered. 

Here we present a statistical parameter extraction test: 
to check for bias in our likelihood code, we construct 500 
realizations of multi-year, multi-frequency sky map data, 
then fit parameters from each realization independently. 
The input maps are transformed to power spectra and 
a likelihood code using the WMAP flight pipeline. The 
simulated inputs, which are also used in ^ are discussed 
in more detail in Appendix lAl 

To determine if the parameter fits are biased, one could 
sample parameters from a ACDM Markov chain for each 
realization individually, then form a weighted average. 
This would be fine if the recovered parameter likelihoods 
were Gaussian, but that is not guaranteed. The opti- 
mal way to combine the likelihoods is to multiply them, 
and then sample from the joint distribution. The prod- 
uct of these likelihoods would represent what we know 
about the universe if we had access to CMB data from 
500 Hubble volumes. Perhaps the most obvious way to 
sample from this distribution is to run a Markov chain. 
However, each of the 500 likelihood functions involves 
an independent main program, set of data inputs, and 
running environment, so this solution is impractical. 
Our approach is to use importance sampling. We want 
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to draw samples from the joint distribution correspond- 
ing to the sum of = 500 log likelihoods. Importance 
sampling draws samples from a covering distribution that 
is close to the desired distribution; it then weights the 
samples by the ratio of probability densities of the two 
distributions t o correct for the difference between the two 
(|Mackavl2003[ ) . This approach allows us to parallelize the 
processing as follows. We generate M = 10,000 samples 
from the covering distribution, compute the model spec- 
tra for each sample and store them, then for each of N 
copies of the likelihood, we separately calculate the log 
likelihood for each sample spectrum. We add the log 
likelihoods, and subtract the log density of the Gaussian 
at that location, to form a weight for each sample. This 
set of weighted samples is effectively the joint likelihood 
of cosmological parameters over the 500 realizations. It 
is not necessary to load N copies of the likelihood code 
into memory, nor establish inter-process communication. 

Determining a useful covering distribution is an iter- 
ative process. We start with a Gaussian model with no 
correlations between parameters, using only one realiza- 
tion, and gradually add realizations (shrinking the region 
of interest), updating the covariance matrix. In the end, 
we use a Gaussian distribution with a covariance ma- 
trix that is a factor of 2 larger than the covariance of 
the 500 combined likelihoods. This makes the width of 
the distribution \/2 times "too large" in each dimension. 
Making the covariance matrix slightly larger than the 
desired distribution allows us to check for the possibility 
of large tails in the distribution of parameters. Because 
the Gaussian sampling distribution has tails which drop 
off exponentially quickly, it could (in principle) fail to 
properly sample the tail of a distribution which fell less 
rapidly. To test this possibility, an array of 2-D scatter 
plots is made with one cosmological parameter on each 
axis, and with the data points color-coded by weight. Vi- 
sual inspection of these plots shows that the weights of 
the sampled distribution are largest in the center of the 
sampling distribution, so that the sampling distribution 
adequately covers the tails of the likelihood. To verify 
that we were not merely looking at the data in a mis- 
leading projection, we also perform a principal compo- 
nent analysis on the sampled points and redisplay them 
plotted on principal component axes. The weights re- 
main highest in the center of the distribution. 

The likelihood for these simulations differs from that 
used in the WMAP 7-year analysis in the following ways. 
Most of these differences are for computational conve- 
nience. Because of disk space and computational time 
limitations, we do not simulate 500 sets of the time or- 
dered data and map reconstruction; we begin with the 
maps. No foregrounds are included, for simplicity. No 
beam error is included, so we do not have to simulate 
small differences in deconvolution. No unresolved point 
source error is included, and we do not introduce point 
sources into the maps. CAMB is run at slightly higher 
accuracy, but this has a negligible effect on parameters. 
The A'sido = 16 pixel likelihood is used for I < 30 temper- 
ature, instead of the Gibbs likelihood, because while the 
Gibbs likelihood runs more rapidly, the pixel likelihood 
data products are much faster to generate. The SZ effect 
is not included in the input spectrum, and so we do not 
attempt to fit it (unlike the 7- year Markov chain analy- 
sis, discussed in M.l^ . While these difi'erences mean that 



the parameter recovery sims do not simulate every part 
of the data analysis, they simulate a substantial portion. 

Since we have 500 times more data in these parameter 
recovery simulations than in the 7- year WMAP data, a 
level of bias well below la is detectable in the simula- 
tions, whereas it is not in the WMAP fiight data. The 
samples from the importance sampling can be fit with 
a 6-dimensional Gaussian. Using this Gaussian to ap- 
proximate the joint likelihood of all 500 realizations, the 
input parameters have = 27 for 6 degrees of freedom, 
which indicates a strong detection of a difference between 
the input and recovered parameters. The assumption of 
Gaussianity is reasonably good here, but the detection 
can also be stated without that assumption. Suppose 
one draws a line (in this 6-dimensional parameter space) 
between the input parameters and the mean of the recov- 
ered parameters, and then projects all the sampled points 
onto that line. Then one can determine how far away the 
input parameters are from the recovered parameters by 
calculating how much weight is on either side of the line 
from the input parameters. In this case, all of the weight 
except for one point is one one side of the line. Because 
that one point is out in the tail of the distribution, it has 
a weight well below average, and much less than one part 
in 10,000 of the weight is on the far side of the input pa- 
rameters. This indicates that the input parameters are 
not consistent with the likelihood of recovered parame- 
ters. Both of these arguments indicate that the input 
parameters are biased. 

However, this analysis shows that the recovered param- 
eters have very little bias compared to their uncertainties. 
The measured level is less than 15% of the 7-year error 
on each parameter, and complete results are given in Ta- 
ble [2] Since the magnitude of the bias is small compared 
to the errors in the 7-year WMAP data, and since it will 
be different for different cosmological models, we do not 
attempt to remove it from the recovered parameters. 

Table 2 shows that we tend to overestimate rig by 0.09 
a. Using a different likelihood code on the 5-year data, 
iRudiord et al.l (|2009( ) found a value of that was 0.6 
a higher than our 5-year result, but with the same pre- 
dicted uncertainty. Since we have demonstrated that our 
likelihood is not biased at this level, we must either con- 
clude that a) there is some bias in the Rudjord et al. 
form of the likelihood, b) there is some undetected resid- 
ual bias in our likelihood, or c) that both are (practi- 
cally) unbiased and that different likelihood approxima- 
tions can lead to parameter estimates that differ by this 
magnitude. To resolve this question, one should evaluate 
both likelihood functions on common data simulations, 
then jointly study the performance of the derived pa- 
rameter ensembles. Our parameter recovery simulations 
show that the WMAP likelihood produces only a small 
bias when averaged over many CMB realizations. This 
does not imply that it is a better approximation to the 
exact likelihood than, e.g., the Rudjord et al. form, thus 
a joint comparison over many data realizations would be 
useful. To get a very rough sense of how large a differ- 
ence one might expect, in case c), above, we consider a 
toy model in which we estimate the variance of N random 
numbers from two partially overlapping subsets of the N 
numbers. In a case where 10% of the total data sample is 
disjoint (i.e., 5% in each subsample is independent) the 
two estimates of the parent variance differ, statistically. 
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Table 2 

Parameter Recovery Bias Test" 



Parameter'' 


Input value'' 


Measured bias"^ 


Bias S/N 


Qa 
ria 

T 


2.2622 

0.11380 

0.72344 

2.4588 

0.9616 

0.08785 


0.0013 ± 0.0025 
-0.00033 ± 0.00025 

0.00167 ± 0.00128 
-0.0160 ± 0.0050 

0.0013 ± 0.0006 
-0.00056 ± 0.00070 


0.02 
-0.06 

0.06 
-0.14 

0.09 
-0.04 



Error bar Accuracy" 



-0.035 
0.074 
0.039 
0.087 
0.020 
0.042 



" Parameter recovery results based on 500 Monte Carlo simulations of 7- year WMAP data 
fit to the 6-parametcr ACDM model. 

' Parameter and its input value in the 500 Monte Carlos realizations. 

Bias measured in the composite likelihood derived from 500 MC realizations, quoted as 
the mean and rms of the marginalized distribution. A positive number indicates that the 
recovered value was higher than the input value. 

Measured bias divided by rms of the marginalized likelihood derived from the WMAP 
data. The bias is less than 15% of the Itr error the 7-year data. 

° Fractional error in reported error bar. We compute the standard deviation of (output - 
input) / (output error), and then subtract 1. For the 150 realizations used (this column 
only), we expect fluctuations of ±0.058 (1-sigma). Our results are compatible with this. 



by 0.3 cr, where a is the rms of each of the subsample 
estimates. Thus, for two hkehhood functions to produce 
parameter estimates that differ by 0.6 a (at 95% confi- 
dence), the two functions must, in effect, be re- weighting 
5% of the data. Given the similar construction of the 
two likehhood functions, this seems unhkely, so further 
study wiU be required to understand this difference. 

We can also use the parameter recovery simulations 
to verify that our error estimates are correct. For each 
of 150 of the data realizations, we use a Markov chain 
to compute the mean and 68% confidence interval for 
each parameter. We then examine the distribution of 
the quantity (output value - input value) /(output error) 
which should have a unit variance. The results are shown 
in the last column of Table [U where we find that the 
errors predicted by the Markov Chain agree with the 
true errors, to within the noise expected from the limited 
number of realizations. 

4. COSMOLOGICAL PARAMETERS FROM WMAP 

In this section we discuss the determination of cos- 
mological parameters using only the 7-year WMAP 
data. The measurements obtained by combining 7- 
year WMAP data with other cosmological data sets are 
presented in iKomatsu et al.l (|2010[ ). Our analysis em- 
ploys the same Monte Carlo Mark ov Chain (MCMC] 
formalism used in p revious analys e s jS pergel et al.ll2003l: 
Verde et al.l [200l iSoergel et al.l [2007. : .Dunklev et alJ 
2009HKomatsu et al.ll2009( ). The MCMC formalism nat- 
urally produces parameter likelihoods that are marginal- 
ized over all other fit parameters in the model. Through- 
out this paper, we quote best-fit values as the mean of 
the marginalized likelihood, unless otherwise stated (e.g., 
upper limits). Lower and upper error limits correspond 
to the 16% and 84% points in the marginalized cumula- 
tive distribution, unless otherwise stated. 

4.1. Six-Parameter ACDM 

The ACDM parameters used are the same as in ^ 
and mentioned in Table [TJ except that A^t, is now 
also sampled. This is a scal e factor for the predicte d 
Sunyaev-Zel'dovich spectrum (jKomatsu fc Selia3l2002D . 
measured at V band, which we add to the TT power spec- 



trum as in iSpergel et al.l (|2007[ ). In the Markov chains, 
this parameter is given a flat prior < Agz < 2, but is 
unconstrained by the WMAP data, so its posterior dis- 
tribution is very fiat over this region. Failing to include 
the SZ effect does not significantly raise the of the fit, 
so only 6 parameters are needed to provide a good fit to 
the WMAP power spectra, and we sample A^z only to 
marginalize over it. 

The ACDM parameters best fit to the 7-year WMAP 
data are given in Table [3l which also lists values derived 
from the 5-year data for comparison. The results are 
consistent, with the 7-year measurements giving smaller 
uncertainties, as expected. The parameters that show 
the greatest improvement are those that most depend on 
the amplitude of the third acoustic peak and the low- 
/ EE polarization: ilbh^, Vtch? , and r, all of which are 
measured about 12% more precisely. The derived late- 
time matter fluctuation amplitude, erg (which depends on 
Oc/i^ and r), is measured 17% more precisely by the new 
data. In ^A.i\ we consider the overall change in allowable 
parameter-space volume offered by the 7-year data. 

As discussed in J}5l this basic ACDM model continues 
to fit the 7-year WMAP data quite well. Indeed, none 
of the additional parameters considered below provide a 
statistically better fit to the 7-year WMAP data, after 
accounting for the fewer degrees of freedom in the fits. 

4.2. Extended Cosmological Models 

In this section we examine the constraints that can 
be placed on augmented ACDM models (and one non-A 
model). In the first group we consider parameters that 
introduce "new physics" : tensor modes, a running spec- 
tral index, isocurvature modes, spatial curvature, and 
non-A dark energy. In the second group, we relax the 
constraints on "standard physics" by allowing the effec- 
tive neutrino number & the primordial helium abundance 
to vary. We also allow the reionization profile to vary. 

4.2.1. Gravitational Waves 

The amplitude of tensor modes, or gravitational waves, 
in the early universe may be written as 



A^(fc) 



27r2 



(9) 
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Table 3 

Six-Parameter ACDM Fit ^ 



Parameter 



7-ycar Fit 



5-ycar Fit 



Fit parameters 





9 9Ka+0-057 

Z.ZOO_g ggg 


2.273 ±0.062 




0.1109 ± 0.0056 


0.1099 ±0.0062 


Ha 


0.734 ± 0.029 


0.742 ±0.030 


is? 


(2.43 ±0.11) X 10"^ 


(2.41 ±0.11) X 10 


Us 


0.963 ± 0.014 


U.yDJ_o 015 


T 


0.088 ± 0.015 


0.087 ±0.017 



Derived parameters 

to 13.75 ± 0.13 Gyr 

Ho 71.0 ± 2.5 km/s/Mpc 

o-g 0.801 ± 0.030 

Hi, 0.0449 ± 0.0028 

fJc 0.222 ± 0.026 

Zeq 31961^!^;^ 

2roion 10.5 ± 1.2 



13.69 ±0.13 Gyr 
71.912 7 km/s/Mpc 
0.796 ±0.036 
0.0441 ± 0.0030 
0.214 ±0.027 
3176li^J 

11.0 ± 1.4 



Models fit to WMAP data only See lKomatsu et al.l 1 12010) 
for additional constraints. 




0.60 0.65 0.70 0.75 0.80 0.85 

Figure 13. Gravitational wave constraints from the 7-ycar WMAP data, expressed in terms of the tcnsor-to-scalar ratio, r. The red 
contours show the 68% and 95% confidence regions for r compared to each of the 6 ACDM parameters using the 7-year data. The blue 
contours are the corresponding 5-year results. We do not detect gravitational waves with the new data; when we marginalize over the 6 
ACDM parameters the 7- year limit is r < 0.36 (95% CL), compared to the 5-ye ar limit of r < 0.43 (95% CL). Tighter limits apply when 
WMAP data are combined with Hq and BAO constraints IIKomatsu et al.ll2010l) . 



where Ph{k) is the power spectrum of tensor perturba- 
tions at wave number k, and t he no rmahzation of Ph{k) 
is as given by iKomatsu et al] ()2009( ). This form is com- 
parable to the curvature perturbation ampHtude, 



~ 27r2 ■ 



(10) 



The dimensionless tensor-to-scalar ratio is defined as 



AUk) 



(11) 



evaluated aX k = 0.002 Mpc~^. In the Markov chain, we 
set a flat prior on r, and require r > 0. 

We do not detect gravitational waves from inflation 
with the 7-year WMAP data, however the upper lim- 
its are 16% lower: r < 0.36 (95% CL) compared to 
r < 0.43 (95% CL). Figure US] shows the the 2-d like- 
lihood contours for r vs. the other ACDM parameters 
using both the 5-year and 7-year WMAP data. This 
shows both the improved upper limit on r and the cor- 
relations with the other measured parameters, especially 
the matter densities and n^. The limits quoted above 
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arise from all of the power spectra measured by WMAP 
with the great est power coming from the shape of the 
TT spectrum. iKomatsu et al.l ()2010[ ) consider the con- 
straints that arise from polarization alone and show that 
the limits improve from r < 1.6 to r < 0.93 using the 
5-year and 7- year data, respectively. 

4.2.2. Scale Dependent Spectral Index 

Some inflation models predict a scale dependence 
or "running" in the (nearly) power-law spectrum of 
scalar perturbations. This is conveniently parameter- 
ized by the logarithmic derivative of the spectral index, 
dris/dhik, which gives ris e to a spectrum of the form 
(|Kosowskv fc Turnei|[T99l 



A2 



ik) 



A2 



n 



(ko) 



\^ 7is{kQ) — l + ^ \n{k/kQ)d7is/d\nk 

(12) 

In the Markov chain, we use a 



with fco = 0.002 MpC 
flat prior on dris/dlnk. 

We do not detect a statistically significant (i.e., >95% 
CL) deviation from a pure power-law spectrum with the 
7-year WMAP data. The allowed range of dus/dlnk is 
both closer to zero and has a smaller confidence range 
using the 7-year data: dus/dlnk = —0.034± 0.026 com- 
pared to dus/dlnk = — 0.037 ± 0.028 from the 5-year 
data. 

If we allow both tensors and running as additional pri- 
mordial degrees of freedom, the data prefer a slight neg- 
ative running, but still at less than 2cr. The joint con- 
straint on all parameters in this model is significantly 
tighter with the 7-year data (see M.3|) . The 7- year con- 
straints on models with additional power spectrum de- 
grees of freedom are given in Table ID 

4.2.3. Isocurvature Modes 

In addition to adiabatic fluctuations, where different 
species fluctuate in phase to produce curvature fluctua- 
tions, it is possible to have an ovcrdcnsity in one species 
compensate for an underdensity in another without pro- 
ducing a curvature. These entropy, or isocurvature per- 
turbations have a measurable effect on the CMB by shift- 
ing the acoustic peaks in the power spectrum. For cold 
dark matter and photons, we define the field 



Pc 



(13) 



(jBean et al.l I2006t IKomatsu et"ani2009D . The relative 
amplitude of its power spectrum is parameterized by a, 



1 — a 



Psjko) 



(14) 



with ko = 0.002 Mpc-\ 

We consider two types of isocurvature modes: those 
which are completely uncorrelated with the curvature 
modes (with amplitude ao), motivated with the axion 
model, and those which are anti-correlated with the the 
curvature modes (with amplitude a_i), motivated with 
the curvaton model. For the latter, we adopt the con- 
vention in whic h anticorrelation incr eases the power at 
low multipoles (jKomatsu et al.|[2009[ ). For both ao and 
a_i, we adopt a fiat prior and require ao > 0, a-i > 0. 



a 



1 1 1 1 


i 1 1 1 1 i 1 1 i 


1 1 1 1 i 1 1 i 1 1 1 1 i 

ui — on 
pIq = yu 
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- 




Ho = 40 " 




\ 


Ho = 30 
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. . 1 . . . 1 . . 
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Figure 14. The range of non-flat models consistent witti tlie 
WMAP 7-year data. The plotted points are drawn from the 
Markov chain (ACDM, with curvature, fit to WMAP data only); 
the color of each point indicates the Hubble constant for that ele- 
ment in the chain, as indicated in the legend. Due to the geometric 
degeneracy, CMB data alone do not strongly constrain closed mod- 
els with Q\ + Qrn > 1, provided a low Hubble constant is tolerated, 
see Table |6] The dashed line parameterizes the geometric degener- 
acy in the 7-year data: Uk = -0.2654 + 0.3697f7A. When WMAP 
data are combined with Ho and BAO data, the re sult strongly fa- 
vors a flat universe: 0.99 < Qtot < 1-01 (95% CL) l IKomatsu et al.l 
120101) . The joints constraints, shown as 68% and 95% blue con- 
tours, provide compelling support for basic ACDM. 

The constraints on both types of isocurvature modes 
are given in Table [5] We do not detect a significant 
contribution from either type of perturbation in the 7- 
year data. The limit on uncorrelated modes improves the 
most with the new data: from ao < 0.16 (95% CL) to 
ao < 0.13 (95% CL) using the 5-year and 7-year data, re- 
spectively. Table [5] also shows that the standard ACDM 
parameters are only weakly affected by the isocurvature 
degrees of freedom. IKomatsu et al.l (|2010D derive anal- 
ogous constraints using a combination of WMAP plus 
other data. They find limits that are roughly a factor of 
two lower than the WMAP-only limits. 

4.2.4. Spatial Curvature 

The basic ACDM model of the universe is flat, with 
rjfc = 1 — r^tot — 0. There is a compelling theoreti- 
cal case for a flat universe in General Relativity, arising 
from the apparent paradox that a flat geometry is dy- 
namically unstable. That is, in order for the curvature 
to be acceptably small today, say < 1, the curva- 
ture in the early universe had to be extraordinarily fine 
tuned. Cosmological inflation achieves this by expanding 
the primordial curvature scale, if any, to super-horizon 
scales today. 

With knowledge of the redshift of matter-radiation 
equality, the acoustic scale can be accurately computed 
for use as a standard ruler at the epoch of recombina- 
tion. The flrst acoustic peak in the CMB then provides 
a means to measure the angular diameter of the acoustic 
scale at the surface of last scattering. If we have indepen- 
dent knowledge of the local distance-redshift relation (the 
Hubble constant, Hq) we can infer the physical distance 
to the last scattering surface, and hence the geometry 
of the universe. If we assume nothing about Hq we are 
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Table 4 

Primordial Power Spectrum Constraints 



Parameter 



ACDM+Tensors 



ACDM+Running ACDM+Tensors+Running 



Fit parameters 



n, 

T 

r 



+0.00073 



0.02313_go(,Q-,2 
U.lUb8_o 0053 
0.757 ±0.031 
(2.28 ±0.15) X 10-9 



0.982 



-0.020 



-0.019 

0.091 ±0.015 



< 0.36 (95% CL) 
dus /dink 
Derived parameters 

to 13.63 ± 0.16 Gyr 

Ho 73.5 ± 3.2 km/s/Mpc 

0-8 0.787 ± 0.033 



n2185+'' °°°^2 

0.688tHS 
(2.42 ±0.11) X lO"'' 
-I p,r,7+0.050 
^■'^^'-0.051 

0.092 ±0.015 



-0.034 ±0.026 

13.87tO-/i^ Gyr 
67.5 ±3.8 km/s/Mpc 
0.818 ± 0.033 



n 09991 +0 00085 
u.uzzi;i_g ooogg 

n 1 1 K7+0.0086 
-0.050 

{2.23t°o:\l) X 10-9 
1.076 ± 0.065 
0.096 ±0.016 

< 0.49 (95% CL) 
-0.048 ± 0.029 



13.79 ±0.18 Gyr 
ig.ll*'" km/s/Mpc 
0.808 ± 0.035 



Models fit to 7-year WMAP data 
straints. 



only. See IKomatsu et al.l II2010I) for additional con- 



Table 5 

Constraints on Isocurvature Modes 



Parameter 



ACDM'' 



ACDM+anti-correlated= ACDM+uncorrelated'* 



Fit parameters 



n n99t;s+0-0005'7 

0.1109 ±0.0056 
0.734 ± 0.029 
(2.43 ±0.11) X 10-9 
0.963 ±0.014 
0.088 ±0.015 



Qa 

n, 

T 

a_i 

Derived parameters 

to 13.75 ±0.13 Gyr 

Ho 71.0 ± 2.5 km/s/Mpc 

as 0.801 ± 0.030 



n n99Q':!+0-00060 

n int;s+0-0057 
o.ioa»_Q 0Q5g 

0.766 ± 0.028 

(2.24 ±0.13) X 10-9 

0.984 ±0.017 

0.088 ±0.015 

< 0.011 (95% CL) 



13.58 ±0.15 Gyr 
74.5^3;^ km/s/Mpc 
0.784-* 



1-1-0.033 
-0.032 



n 0231 5+°-9oo''^ 

0.758 ±0.030 
(2.38 ±0.11) X 10-9 
0.982 ±0.020 
0.089 ±0.015 



< 0.13 (95% CL) 

13.62 ±0.16 Gyr 
73.6 ±3.2 km/s/Mpc 
0.785 ± 0.032 



^ Models fit to 7- year WMAP data only. See IKomatsu et al.l l|201CII ') for additional con- 
straints. 

^ Repeated from Table [3] for comparison. 

Adds curvaton-type isocurvature perturbations IjKomatsu et al.ll2010l '). 
Adds axion-type isocurvature perturbations IIKomatsu et al.ll2010l) . 



left with a geometric degeneracy which is iUustrated in 
Figure [Til 

Assuming ACDM dynamics, WMAP data alone pro- 
vide a remarkably simple constraint on the geometry 
and matter-energy content in the universe. The geo- 
metric degeneracy in the 7-year data is well described 
by Qk = -0.2654 + 0.3697f^A (the dashed line in Fig- 
ure [T4|). We have placed a flat prior on ftk, and we now 
also constrain fl\ > 0. The Figure also quantitatively 
illustrates how knowledge of the Hubble constant fixes 
the geometry, + Qm, and vice- versa. The points in 
the plot arc culled from the Markov Chain that samples 
this model and their color is coded by the value of the 
Hubble constant for that sample. As one moves down 
the degeneracy line, the Hubble constant must decrease 
for the model to remain consistent with the geometry 
imposed by the CMB. For a flat universe the 7- year data 
give Ho = 71.0 ± 2.5 km/s/Mpc (Tabled, in excellent 
agreement with more traditional measurements of the 



Hubble constant, e.g. lRiess et al.l (j2009t ). 

If we allow curvature as a parameter, the 7-year 
WMAP data improve on the 5-year constraint by 11%, to 
Qk ~ — 0.080lQ Qg3. While this result is consistent with 
a flat universe, the preferred model is slightly closed and 
has a relatively low Hubble constant due to the geomet- 
ric degeneracy, Hq ~ 531^5 km/s/Mpc. Therefore, if we 
impose local distance scale measurements in the form of 
Hq and BAO data, the limits on curvature tighten signif- 
icantly to rik = -0.00231°;°°^^ (IKomatsu et al.l[2010l) . 

4.2.5. Non-A Dark Energy 

Dark energy is believed to be driving the present-day 
acceleration of the universe. Current measurements are 
consistent with the dark energy being a cosmological con- 
stant, or vacuum energy. If it is not a cosmological con- 
stant, then its physical density may change with the ex- 
pansion of the universe. This, in turn, would affect the 
expansion history and the rate of large-scale structure 
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Table 6 

Constraints on Curvature and Dark Energy ' 



Parameter 



ACDM'' 



OACDM': 



vCDW^ 



Fit parameters 



n n99c;s+0 00057 



Ha 



0.1109 ±0.0056 
0.734 ±0.029 
(2.43 ±0.11) X 10-9 
ris 0.963 ±0.014 

T 0.088 ±0.015 

w 

Derived parameters 

to 13.75 ± 0.13 Gyr 

Ho 71.0 ± 2.5 km/s/Mpc 

as 0.801 ± 0.030 



0.02229 



+0.00058 
-0.00057 



0.1117"* 



7+0.0053 
-0.0055 

< 0.77 (95% CL) 
(2.48 ±0.11) X lO-** 
0.955 ± 0.014 
0.086 ± 0.015 



-0.080 



+0.071 



53^ 



15.9+^:1^ Gyr 

km/s/Mpc 
762+° ''*'' 



a+0. 00063 
-0.00062 

0.1112 ±0.0058 



0.02258^ 



0.741 



+0.095 



-0.099 

(2.43 ±0.12) X 10" 
0.964 ±0.015 

;i + 0.016 

-0.015 



0.088"* 



-1.12 



+0.42 
0.43 



13.75tO-J9 Gyr 
75^11 km/s/Mpc 
0.83 ±0.14 



Models fit to 7- year WMAP data only. See IKomatsu eFall H2010I 1 for additional 
constraints. 

^ Repeated from Table [3] for comparison. 

Adds spatial curvature as a parameter, writh w = —I. Q\ > is imposed as a 
prior. 

Adds dark energy equation of state as a parameter, with fi^ = 0. w > —2.5 is 
imposed as a prior. 



growth in the universe. The evolution of its physical 
density is governed by its equation of state w = p/p 
where p is the pressure of the dark energy and p its den- 
sity. The cosmological constant has an equation of state 
w = — I. It would be tremendously important if observa- 
tions could determine that w ^ —1 since that would rule 
out the prime candidate for the dark energy and provide 
important new clues about physics. 

Since the CMB primarily probes the high redshift uni- 
verse {z ^ 1000), and the effects of dark energy only 
start to dominate at relatively low redshift (z ~ 2), the 
CMB is not especially sensitive to subtle properties of 
the dark energy. Nonetheless, meaningful constraints on 
the equation of state can be inferred from the 7-year 
WMAP data. If we assume the universe is flat but let 
w be a parameter in the Friedmann equation (with a 
flat prior on w, —2.5 < w < 0, and w' = 0), we ob- 
tain the constraints given in Table |6l In particular, the 
7-year data give w = —1.12 '^'^'^^, which is consistent 
with a cosmological constant. IKomatsu et al.l ()2010[ ) in- 
vestigate the constraints imposed when 7-year WMAP 
data are combined with other observations. With BAO 
and Hubble constant measurements added, they find 
w = —1.10 ±0.14, which provides compelling limits on w 
without using Type la supcrnovae data. When SNe data 
are included, the result becomes w = —0.980 ± 0.053, 
but the quoted error does not include systematic errors 
in the supernovae, which are comparable to the statistical 
errors. Accounting for this would produce a final uncer- 
tainty that is roughly half the size of the error without 
SNe data. 

If we relax the assumption tha t w' = and/or S7fc = 
0, the constraints on w weaken (|Komatsu et al.|[2010[ ). 
This points to the need for more accurate and precise 
measurements of the expansion history and growth rate 
of structure if we are to gain further clues about dark 
energy from cosmology. 



4.2.6. Neutrinos 

Neutrinos affect the CMB spectrum in a variety of 
ways; one is by providing relativistic degrees of freedom 
to the plasma prior to recombination. Since neutrinos, or 
other relativistic species, are not coupled to the photon- 
baryon fluid, they free-stream out of over-densities and 
damp the acoustic oscillations prior to recombination. 
This action suppresses the peaks in the angular power 
spectrum somewhat; the amplitude of the effect depends 
on the effective number of relativistic degrees of freedom. 
Using the 7-year WMAP data we place a 95% CL lower 
limit of iVcff > 2.7 (95% CL) on the number of relativis- 
tic degrees of freedom, for a flat prior < Ngg < 10. 
(The standard model has Ncs = 3.04.) This new limit is 
17% higher than the 5-year limit of iVcff > 2.3 (95% CL) 
- due to the improved third peak measurement - and is 
now quite close to the standard model value. 

The mean energy of a relativistic neutrino at the epoch 
of recombination is (E) = 0.58 eV. In order for the CMB 
power spectrum to be sensitive to a non-zero neutrino 
mass, at least one species of n eutrino must have a m ass 
in excess of this mean energy (jKomatsu et al.l [20091 ). If 
one assumes that there are N^s = 3.04 neutrino species 
with degenerate mass eigenstates, this would suggest 
that the lowest total mass that could be detected with 
CMB data is T^.jrii, ^ 1.8 eV. Using a refined argument, 
llchikawa etHTlpOOSl ) argue that one could reach ~1.5 
eV. When we add ^ as a parameter to the ACDM 
model (and a flat prior on the physical neutrino density, 
fl,yh^, constrained by fl^h^ > 0), we obtain the fit given 
in Table [71 specifically < 1.3 eV (95% CL), which 

is unchanged from the 5-year result and is slightly below 
the basic limits just presented. Note that these results 
come from the WMAP data alone. Tighter limits may 
be obtained by combining CMB data wit h measurements 
of str ucture formation, as discussed in IKomatsu et al.l 
(|2010f l and references therein. 
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Table 7 

Constraints on Neutrino Properties ^ 



Parameter 



ACDM'^ 



ACDM+JVeff'^ 



ACDM+E 



Fit parameters 
Ha 



0.1109 ±0.0056 



A2 



0.734 ±0.029 
(2.43 ±0.11) X 10~^ 
lis 0.963 ± 0.014 

T 0.088 ±0.015 

E "^-^ 

Derived parameters 

to 13.75 ±0.13 Gyr 

Ho 71.0 ± 2.5 km/s/Mpc 

(T8 0.801 ± 0.030 



0.02260 



+0.00055 
-0.00054 
,+0.041 



0-162lo.o38 

r, 701+0.029 

(2.39 ±0.11) X 10- 
002+0-022 

0.088lHlt 
> 2.7 (95% CL) 



11.9: 

s+8.9 



Gyr 



3+0.00061 
-0.00062 

0.1122 ±0.0055 



0.02219^ 



660+°"^2 
(2.50 ±0.13) X 10-0 
0.953 ±0.017 
0.086 ±0.014 

< 1.3 eV (95% CL) 



82.6I0 7 km/s/Mpc 65.014 5 km/s/Mpc 



-8.7 

0.903 ±0.077 



14.09 
+4.4 



+0.25 



Gyr 



4.5 

0.685 



+0.079 



Models fit to 7- year WMAP data only. See IKomatsu eFall II2010I ') for additional 
constraints. 

^ Repeated from Tabic |3] for comparison. 

Allows effective number of relativistic species to vary (cf. A^cff=3.04). N^g < 10 
is imposed as a prior. 

Adds neutrino mass, J^m^, as a parameter, assuming A^(,ff=3.04 and degenerate 
mass eigenstates. 



4.2.7. Width of Reiomzation 

Effecti ve with the March 2008 version of the code 
CAMB (iLewisI 120081 ). a new parameter has been added 
which allows users to vary the reionization profile while 
holding the total optical depth fixed. The basic profile 
is a smooth ramp in redshift space and the parameter, 
Az, changes the slope of the ramp about its midpoint in 
such a way as to preserve total optical depth. 

We have added as a parameter to the basic ACDM 
model, with a flat prior in the range 0.5 < < 15, and 
present the results in Table |S1 

4.2.8. Primordial Helium Abundance 

Helium is thought to be synthesized in the early 
universe via Big Bang Nucleosynthesis (BBN). Given 
the WMAP measurement of the baryon-to-photon ra- 
tio, 77, the BBN -predicted yield for helium is FHe=0.249 
()Steigmanll2007[ ). To date, the best technique for measur- 
ing the primordial abundance has been to observe stars in 
HII regions: in these systems, the helium abundance as a 
function of metallieity can be observed and the relation 
can be regressed to zero metallieity, w hich is presumed 
to g iv e the primordia l abu n dance ([Gruenwald et al.l 
200a Hiotov fc ThuanI [2003: lOlive fc SkillmanI 12005 
Fukugita fc Kawasaki"2006!: iPeimbert et al.ll2(307f ) as re- 
viewed by [Steigman (2007). 

Primordial helium affects the time profile of recombina- 
tion which, in turn, affects the CMB angular power spec- 
trum, especially the third acoustic peak. With WMAP's 
improved measurement of this peak, it is now possible to 
let Yho be a fitted parameter in the ACDM model. We 
use a fiat prior with 0.01 < Yho < 0.8. We present the 
results of this fit in Table |8] and call out the helium abun- 
dance specifically: Yhc = 0.28j^Q |^5. This result is consis- 
tent with the BBN prediction and sugges ts the existence 
of pre -stellar helium, at the ~ 2a level. IKomatsu et al.l 
()2010D consider the constraints that can be applied when 
higher-resolution CMB data are included in the fit. They 



find that the combined CMB data produce, for the first 
time, evidence for pre-stellar helium at > 3cr. 

4.3. Volume Change 

The tables presented in this section quote parameter 
uncertainties for marginalized 1-dimensional likelihood 
profiles. With the 7- year and 5-year results side by 
side, one could infer the improvement in precision for 
any given parameter in any given model fit, and we have 
called out examples in the text. But it is difficult to mea- 
sure the overall improvement in the model fits from this 
presentation. A better measure is given by comparing 
the allowable volume in A'^-dimensional parameter space 
for each of the models. 

As a proxy for allowable volume, we compute the 
square root of the determinant of the parameter co- 
variance matrix for each model, using data from the 
Markov chains. For example, with the 6-parameter 
ACDM model, we compute the 6x6 eovarianee matrix 
directly from the chain samples. The allowable volume 
ratio is defined as the ratio of the square root of the 
5-year determinant to the corresponding 7-year value. 
Models with more fit parameters have more rows and 
columns in their eovarianee matrix. While this proxy 
is only proportional to the volume if the parameter dis- 
tributions are Gaussian, the error in this approximation 
will be similar for both the 5-year and 7-year data sets, 
so a comparison is still valid. 

Table [9] gives the change in allowable parameter-space 
volume as a ratio of the 5-year to 7-year value. These 
results are based on fits to WMAP data only. Overall, 
the 6-parameter ACDM model is measured a factor of 1.5 
more precisely with the 7- year data while the model with 
2 additional parameters, tensors plus a running spectral 
index, is measured a factor of 3 times more precisely. 
Models with one additional parameter typically improve 
by factors of 1.8 to 1.9. 
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Table 8 

Tests of Standard Model Assumptions^ 



Parameter ACDM*' ACDM+A^'^ ACDM+ync'^ 

Fit parameters 

Uch^ 0.1109 ±0.0056 

f^A 0.734 ±0.029 

A^, (2.43 ± 0.11) X 10-3 

Us 0.963 ±0.014 

r 0.088 ±0.015 

Derived parameters 

to 13.75 ± 0.13 Gyr 

Ho 71.0 ± 2.5 km/s/Mpc 

as 0.801 ± 0.030 

Models fit to 7- year WMAP data only. See IKomatsu eFall 1)20101 ) for additional 
constraints. 

^ Repeated from Table |3] for comparison. 

Allows width of reionization to vary. The constraints on arc limited by the prior 
0.5 < A2 < 15 and so are not listed in the table. 

Allows primordial helium mass fraction to vary (cf. Yho=0.24). 



Table 9 

Parameter-Space Volume Reduction^ 



Model 


Dimension 


Ratiol" 


ACDM (Table [Sj 


6 


1.5 


ACDM+tensors (Tabled 


7 


1.9 


ACDM+running (Tableg} 


7 


1.7 


ACDM+tensors+running (Table Ul 


8 


3.0 


ACDM+anticorrclated isocurvature (Table [St 


7 


1.9 


ACDM+uncorrelatod isocurvature (Table [5]l 


7 


1.9 


ACDM+massive neutrinos (Table |3|l 


7 


1.8 


oACDM (Tabic Hi 


7 


1.8 


uiCDM (Tabic |6} 


7 


1.5 



^ The relative change in allowable parameter-space volume when mod- 
els are fit to the 7-year WMAP data in place of the 5-year data. For 
a given model, the allowable volume is defined as the square root of 
the determinant of the parameter covariance matrix, as obtained from 
the Markov chains. The basic set of 6 parameters compared in the 
first row are {Cftfc'^, Slc'i^, f^Ai "-s, t}. Additional parameters are 
as noted in the first column. 

^ The ratio of allowable parameter-space volume: 5-year over 7-ycar, 
when fit to WMAP data only. 

5. GOODNESS OF FIT 

Given a best-fit model from the MCMC analysis, we 
can ask how well the model fits the data. Given that 
the likelihood function is non-Gaussian, answering the 
question is not as straightforward as testing the per 
degree of freedom of the best-fit model. Instead we re- 
sort to Monte Carlo simulations and compare the abso- 
lute likelihood obtained from fitting the flight data to an 
ensemble of simulated values. 

For testing goodness of fit, we generate 500 realizations 
of the 7-year sky map data that include the CMB signal 
and instrument noise. These are the same realizations as 
were used for parameter recovery in ^ and are discussed 
in more detail in Appendix [Al 

For each realization of a 7-year data set, we constructed 
the likelihood function appropriate to those data. This 
required forming the high-Z MASTER spectra and the 
low-resolution sky maps used in the code. (For this study 
we did not employ Gibbs sampling for the low-l TT like- 



0.02244 ± 0.00055 0.02253l° o°g5g 

7+0.0054 n 1 1 qn+O-OO^g 



0.728 ±0.028 



-1+0.031 



(2.46 ± 0.12) X 10-9 (2.43 ± 0.11) X 10" 
0.958lHi^ 0.969lHl^ 
0.087 ±0.015 0.0881q;oj5 



0.28 



+0.14 
-0.15 



U.ntoll Gyr 13.69 ± 0.16 Gyr 

70.5 ± 2.4 km/s/Mpc 70.9 ± 2.5 km/s/Mpc 
0.802 ±0.030 0.820loo54 



lihood, rather we used a direct pixel-space code that was 
computationally slower than the Blackwell-Rao estimate 
per likelihood evaluation, but it required less setup over- 
head per data realization.) For goodness of fit testing, 
we evaluated the likelihood of the input ACDM model 
for each data realization. 

Due to its hybrid nature, the likelihood produces sev- 
eral components that need to be combined to obtain the 
full likelihood. The components of most interest to good- 
ness of fit testing are the high-Z TT and TE portions, 
which cover the bulk of the multipole range and are the 
most straightforward to interpret. Recall the the high- 
/ TT component contains both a Gaussia n and a log- 
norma l contribution, as per equation 11 of iVerde et al.l 
(|200l 

1 2 

ln£ = -In^Gauss + 2 In^LN- (15) 

Here the first term can be compared to y^, as per equa- 
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tion 6 of IVerde et all (|2003l) 

In /:Gauss (X -i ^(C, - Cl'')Qu' (CV - 4'^). (16) 
w 

where Q is the inverse covariance matrix of the observed 
power spectrum C, and C**^ is the model spectrum. The 
high-Z portion of the TE hkelihood includes only the 
Gaussian component, which is a good approximation 
given the lower signal-to-noise ratio of the TE data. In 
the following, we report on the distribution of — 21n£ 
which we call the effective x^- 

We compare the distribution of effective values for 
the high-Z TT portion of the likelihood, which contains 
1170 multipoles from I = 31-1200, to a distribution 
with 1170 degrees of freedom. The agreement between 
the two distributions is good. We tentatively attribute a 
small shape difference between them to the non-Gaussian 
component in the likelihood. The effective for the 7- 
year flight TT spectrum is 1227 for 1170 degrees of free- 
dom, after marginalizing over point sources and the SZ 
spectrum (which arc not in the simulations). According 
to the Monte Carlo distribution, 48/500 of the realiza- 
tions had a higher effective x^, indicating that the flight 
data are reasonably well fit by the ACDM model spec- 
trum. 

We perform a similar comparison for the high-/ TE 
data, which covers the multipole range I = 24-800. One 
point of note is that we have adjusted our empirical cal- 
ibration of /sky,TE as a result of these simulations: the 
new value is 1.011 times larger than we used in the 5- 
year analysis, which is equivalent to over-estimating the 
5-year TE errors by 1.1%. The re-calibrated Monte Carlo 
distribution tracks the pure x^ distribution, consistent 
with the high-/ TE likelihood being Gaussian. The ef- 
fective x^ for the 7-year flight TE spectrum is 807 for 
777 degrees of freedom (again, after marginalization over 
point sources and the SZ spectrum). According to the 
Monte Carlo distribution, 113/500 of the simulations had 
a higher effective x^, which could easily happen by ran- 
dom chance. This indicates that the ACDM theory yields 
a TE spectrum that fits the data well. 

6. CONCLUSIONS 

We present the angular power spectra derived from the 
7-year WMAP sky maps and discuss the cosmological 
conclusions that can be inferred from WMAP data alone. 

With the 7-year data, the temperature (TT) spec- 
trum measurement is now limited by cosmic variance 
for multipoles / < 548, and the signal-to-noise ratio 
per multipole exceeds unity for / < 919. In a band- 
power of width Al = 10 the signal-to-noise ratio exceeds 
unity to / = 1060. The third acoustic peak in the TT 
spectrum is now well measured by WMAP. In the con- 
text of a flat ACDM model, this improvement allows us 
to place tighter constraints on the matter density from 
WMAP data alone, Qmh"^ = OASSAtomlt and on the 



epoch of matter-radiation equality, Zgq = 3196j;i[33. The 
temperature-polarization (TE) spectrum is detected in 
the 7-year data with a significance of 20(T, compared to 
13(7 with the 5-year data. We now detect the second 
dip in the TE spectrum near / ~ 450 with high confi- 
dence. The TB and EB spectra remain consistent with 
zero in the 7- year data. This demonstrates low system- 
atic errors in the data and is used to place 33% tighter 
limits on the rotation of linear polarization due to parity- 
yiolating effects : Aa = -1.1° ± 1.4°(stat.) ± 1.5°(sys.) 
(jKomatsu et al.ll2010[ ). The low-/ EE spectrum, a mea- 
sure of the optical depth due to reionization, is de- 
tected at 5.5(7 significance when averaged over / = 2- 
7: /(/ + l)Cf^/{2Tr) = 0.074t°:°^4 (68% CL). 

The high-/ EE spectrum in the range 24 < / < 800 
is detected at over 8 a. The BB spectrum, an im- 
portant probe of gravitational waves from inflation, re- 
mains consistent with zero; when averaged over / = 2-7, 
/(/ + l)Cf ^/(27r) < 0.055 /iK^ (95% CL). The upper 
limit on tensor modes from polarization data alone is a 
factor of 2 lower with the 7-year data than it was using 
the 5-year data (jKomatsu et al.ll2010l ). 

The data remain consistent with the simple ACDM 
model. The best-fit ACDM parameter values are given 
in Table [3] the TT spectrum from this fit has an ef- 
fective x^ of 1227 for 1170 degrees of freedom, with a 
probability to exceed of 9.6%. The allowable volume in 
the 6-dimensional space of ACDM parameters has been 
reduced by a factor of 1.5 relative to the 5-year vol- 
ume. Most models with one additional parameter beyond 
ACDM see volume reduction factors of 1.8-1.9, while the 
ACDM model that allows for tensor modes and a run- 
ning scalar spectral index has a factor of 3 lower volume 
when fit to the 7-year data. We test the parameter re- 
covery process for bias and find that the scalar spectral 
index, Ug, is biased high, but only by 0.09a, while the 
remaining parameters are biased by < 0.15a. 

The improvement in the third peak measurement leads 
to tighter lower limits from WMAP on the number of rel- 
ativistic degrees of freedom (e.g., neutrinos) in the early 
universe: iVeff > 2.7 (95% CL). Also, using WMAP data 
alone, the primordial helium mass fraction is found to be 
Ykc = 0-28tQ;i[g, and with data from higher-resolution 
CMB experiments included, IKomatsu et al.l (120101 ) es- 
tablish the existence of pre-stellar helium at > 3(7. 
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APPENDIX 

A. PARAMETER RECOVERY SIMULATIONS 



This appendix describes the configuration of the set of 500 simulations that was used for checking parameter recovery 
and the x^ values from the likelihood. 
The mask is the KQ85y7 mask, which lets through 78.3% of the sky. 



WMAP SEVEN-YEAR OBSERVATIONS: POWER SPECTRA AND WMAP-DERIVED PARAMETERS 21 



For instrument noise, we employ a 2-step process in which we generate uncorrelated noise at high resolution and 
combine it with low resolution correlated noise. The noise was constructed from the 7-year iVsido = 512 and 1024 A^obs 
maps, t he 7-year A^s i rie = 16 covariance matrices, the 7-year ctq values, and the 7-year synchrotron cleaning factors, 
given in iGold et al.l ()2010t ). For each year and nine DAs (Kal-W4), a high resolution uncorrelated noise map was 
made with I, Q, U components at A^sido = 512 and and 1024, including the QU correlations within each pixel. The 
noise maps are generated on a single-year, single-DA basis so wc can mimic the construction of the flight spectra and 
likelihood function. A correlated noise map was made at Ngido = 16. The ctq values for the polarization portion of 
these maps were increased by a factor of 1/(1 — oi), where oi is the fraction of the K band map which has been 
removed to avoid synchrotron contamination ([Gold et al.ll2010[) . This accounts for the increased noise due to the 
synchrotron template subtraction. We combined the low resolution (A^sido = 16) and high resolution (A^sidc = 512 or 
1024) maps by subtracting off the mean of the high resolution noise within each low resolution pixel, and then added 
the low resolution noise to all high resolution pixels within that low resolution pixel. This process provides a high 
resolution noise realization that has the proper low resolution noise correlations when it is binned. Note that the 
A'sidG = 512 and 1024 maps have different noise realizations, but this will have no effect on the resultant likelihoods, 
because the A^sido — 512 maps are used for polarization and the A'sidc = 1024 maps are used for temperature, following 
the procedure used in the standard WMAP pipeline. 

The CMB signal is assumed to be Gaussian, and random-phase, and so its statistical properties are completely 
defined by a power spectrum. The parameter recovery simulations all use the same power spectrum, which was 
derived from the best fit to a 5-year ACDM Markov chain, with WMAP, Baryon Acoustic Oscillation, and Supernova 
data. The parameters used are rj^/i^ = 0.0022622, Qch'^ = 0.1138, Hq = 70.234 kms-iMpc"\ A'^ = 2.4588 x lO'^, 
Us = 0.9616, T = 0.08785. The gravitational lensing signal is treated as Gaussian, and the effects on the temperature 
and polarization power spectra are included. However, the BB spectrum has been zeroed, for consistency with a 
map-making simulation done previously. Zeroing the BB spectrum will have no effect on parameter recovery, since 
our simple ACDM model has a tensor to scalar ratio of 0, and therefore an undetectable BB spectrum. Each CMB 
realization has the same theoretical power spectrum, but different cosmic variance. The different DAs in a given 
realization all see the same CMB sky, but with different smoothing, and different noise. The smoothing used for each 
DA is a circular beam response based on the appropriate 7-year beam transfer function ijJarosik et al.|[2bl0() . 

These simulations do not include foregrounds, the Sunyaev-Zel'dovich effect, unresolved point sources, or beam 
uncertainty. 

The analysis of the parameter recovery simulations uses the same parameters as went into the maps, in the case of 
A'obs, beam profiles, and usage of gravitational lensing. However, for the ctq and synchrotron cleaning coefficients, the 
simulations were produced with the 7-year values, and analyzed with the previous 5-year values. Note that this is a 
small difference, because the 7-year values are not very different from the 5-year values. 

The parameter recovery importance sampling used the September 2008 version of CAMB. and so docs not include 
the bug fix for the proton mass error. This has no effect on the simple ACDM model explored here. 

B. RECFAST AND WMAP LIKELIHOOD UPDATES 

As mentioned in ^2.6\ since the original version of this paper, there have been two updates to the parameter 
estimation code: an improvement in RECFAST from version 1.4.2 to version 1.5, which includes an improved model 
of the hydrogen atom; and a small bug fix in the WMAP likelihood code, changing its version from 4.0 to 4.1. 

We have rerun the ACDM parameter fits using the version 4.1 likelihood, and find parameter changes of order 0.1 
a. The effect on parameters of updating to the new version of RECFAST in CAMB is also of order 0.1 a. The largest 
combined changes arc an increase in erg of 0.35 tr, and an increase in the spectral index of 0.26 a; Table fTOl provides 
more detailed information on several original and updated parameter values. 
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